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Stochastic microlensing is a central tool in probing dark matter on galactic scales. From first prin- 



o 

'— ■ ciples, we initiate the development of a mathematical theory of stochastic microlensing. Beginning 

with the random time delay function and associated lensing map, we determine exact expressions for 



the mean and variance of these transformations. In addition, we derive the probability density func- 
tion (p.d.f.) of a random point-mass potential, which form the constituent of a stochastic microlens 
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^sj potential. We characterize the exact p.d.f. of a normalized random time delay function at the origin 

p' 
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large number of point masses if the normalized time delay function was at a general point of the 

o 



showing that it is a shifted gamma distribution, which also holds at leading order in the limit of a 



lens plane. For the large number of point masses limit, we also prove that the asymptotic p.d.f. of 
the random lensing map under a specified scaling converges to a bivariate normal distribution. We 



show analytically that the p.d.f. of the random scaled lensing map at leading order depends on the 
magnitude of the scaled bending angle due purely to point masses as well as demonstrate explicitly 
how this radial symmetry is broken at the next order. Interestingly, we found at leading order a 
formula linking the expectation and variance of the normalized random time delay function to the 
first Betti number of its domain. We also determine an asymptotic p.d.f. for the random bending 
angle vector and find an integral expression for the probability of a lens plane point being near a 
fixed point. Lastly, we show explicitly how the results are affected by location in the lens plane. 
The results of this paper are relevant to the theory of random fields and provide a platform for 
further generalizations as well as analytical limits for checking astrophysical studies of stochastic 



microlensing. 

PACS numbers: 02.30.Nw, 98.62.Sb, 02.50.-r 

Keywords: gravitational lensing, random mappings, probability, asymptotics 



I. INTRODUCTION 



Stochastic microlensing is the study of the deflection of light by a random collection of stars. In recent years, this 
subject has become a central tool in understanding the nature and distribution of dark matter on galactic scales (e.g., 
Schechter and Wambsganss [1]). This has raised the need for a rigorous framework for stochastic microlensing. 

Research on stochastic microlensing has been mainly numerical, or at least semi-analytical, in flavor. Work in the 
early 1980's included an analytical study of certain stochastic effects of microlensing by a single random lens (e.g., 
Vietri and Ostriker [2||Y. This was later extended to the case of multiple lenses using numerical simulations (e.g., Kayser, 
Refsdal, and Stabell 3|; Paczynski [4J). While successful in furthering our picture of stochastic microlensing, these 
simulations often cannot be checked analytically and are time-consuming. Some of the first analytical counterparts of 
the numerics were focused on the study of random surface brightnesses, bending angles, magnifications, and redshifts 
(e.g., Deguchi and Watson 5|, |fj; Katz, Balbus, and Paczynski 7(; Schneider [§J; Turner, Ostriker, and Gott 9]; 



Vietri 



10(). This was followed by numerical and semi-rigorous work characterizing the magnification probabilities 



12 



and the expected number of microimages in a variety of lensing situations (e.g., Rauch et al. [Ill ]; Wambsganss 
Wambsganss, Witt, and Schneider 13j). Recently, similar techniques were employed to study the distributions of 
microimages, magnifications, and time delay differences between images, with the last two applied directly to probe 
dark matter using flux-ratio anomalies (e.g., Granot, Schechter, and Wambsganss 



14J : Keeton, and Moustakas 



161 ] ; Tuntsov et al. 



m 



171]). In a more general 



Schechter, and Wambsganss [l|; Schechter, Wambsganss, and Lewis 

context, Berry and collaborators also explored the statistical aspects of caustics in optical systems for the case of 

Gaussian random fields (see 18 1 and references therein). 

This series aims to develop a mathematical theory of stochastic microlensing, starting from first principles. This 
fits in the general theory of random fields, namely, the study of random functions and their critical point structure 
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(for more on random fields, see Adler and Taylor 19] and references therein). We derive new stochastic results, for 
a general point in the lens plane, about random time delay functions and lensing maps, which are transformations 
forming the core of microlensing theory. We determine the exact (as opposed to asymptotic) expectation and variance 
of both the lensing map and the point mass potential. Using the expected value of the components of the lensing map, 
we deduce, for a fixed image in the lens plane, the corresponding average position of the point it is mapped to in the 
light source plane. Furthermore, since point mass potentials form the density components of the time delay function, 
we use the exact formulas of their first and second moments to derive the asymptotic expectation and variance of 
the normalized time delay function in the large number of stars limit. We also reveal an interesting topological link, 
namely, at leading order we show that the expectation and variance of the normalized random microlensing time delay 
function is related to the first Betti number of the time delay functions's domain. 

Continuing our study of the normalized random time delay function, we establish that at the origin its exact p.d.f. 
is a shifted gamma distribution, and that this is asymptotically the case for an arbitrary point within the lens plane 
in the large number of stars limit. In the same limit, we prove that the p.d.f of a scaled lensing map converges to 
a bivariate Gaussian, with once again consistent asymptotic mean and variance. This analytical work shows that 
the p.d.f. of the random scaled lensing map at leading order term depends on the square of magnitude of the scaled 
bending angle due purely to point masses. It also illustrate how this radial symmetry is broken at the next order. 
Similarly, we note that the p.d.f. of the random time delay function at the origin depends only on the potential time 
delay due purely to point masses, a result also holding at leading order for an arbitrary point in the lens plane. 

Taking our analysis further, we compute the probability that the image of a given point under the lensing map will 
be in a given disk in light the source plane, which provides information about the light source position. We deduce 
from these results the expectation, variance, and asymptotic p.d.f. of the bending angle vector at an arbitrary position 
in the lens plane. We also use these results to compute the probability that a point in the lens plane will be mapped 
in a certain neighborhood of "itself" in the light source plane. This study relates to the theory of fixed points, which 



was studied earlier by Petters and Wicklin [20| in the case of a deterministic lens system. 

Overall, these first steps lay the foundation in terms of concepts, results, and tools for future analytical work on 
random time delay functions and lensing mappings, including lensing observables like magnification, arising from more 
complex stochastic mass distributions. The work can also help as a guide on tractable aspects of the general theory 
of random mappings. 

We state the basic random microlensing framework in Section [TTJ In Section [ill! we present the results of our study 



of the random time delay function, obtaining its expectation, variance, exact p.d.f. at the origin, and asymptotic p.d.f. 
at an arbitrary lens plane point. A similar presentation of results is done for the random lensing map in Section [IVI 
In Section [Vj we apply our results about random lensing maps to random bending angle vectors and fixed points. 
The detailed proofs of these results are given in the Appendix. 

Finally, the second paper in this program will explore the random shear and expected number of microimages 
produced by stochastic microlensing. 

II. BASIC FRAMEWORK 

Consider a microlensing situation where the light source is point-like and located at y in R 2 (the light source plane) 
and the gravitational lens is a collection of g point masses (stars) that are randomly distributed over a region of the 
lens plane. This lensing scenario is modeled with the following gravitational lens potential, where all quantities are 
in dimensionless form (see 21J, pp. 79, 104]): 

V- fl W = yN 2 -|(x?-xI)+^m i log| a; -^|, (1) 

i=i 

with x = {x\,X2) in the lens plane L = {£i, . . . ,£ ff }. 

In a standard stochastic microlensing situation, the light source position y is assumed to be a uniform random 
vector in R 2 and the star positions £j = (Uj, Vj) are taken to be independent, identically distributed (i.i.d.) uniform 
random vectors in a region of the lens plane. The continuous matter density k c , external shear 7, masses rrij, and 
number of stars g are assumed fixed, unless stated to the contrary. The random lens potential {!]) is widely used in 
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studies of stochastic lensing due to stars (see 
Remark: 

• The constants rrij are scalar multiples of the actual stars' physical masses, but we shall continue to call them 
masses for simplicity (see [2l|, p. 102]). 

• A natural next step is to consider a lensing scenario where the masses are distributed isothermally in an elliptical 
galaxy. This is a work in progress. 

Gravitational microlensing of a light source at position y in R 2 causes a change in a light ray's arrival time that is 



captured analytically by the dimensionless random time delay function given as follows (see 

2 i »\ 2 ' ' 2 



2l|,p. 81): 



T g,v( x ) = t:\x- y\ 2 - -rfM 2 + 77O1 ~ x%) -^mj log\x - £j 



(2) 



21], p. 



i=i 



Fermat's Principle of stationary time states (see |21| . p. 67): Light rays from the source at y to the observer are 
characterized by critical points of the time delay function T v (x), that is, the solutions x in the lens plane L, where 
L = M 2 — {£i, . . . , £ ff }, of the equation 



21], p. 81): 



VT y (x) = 0, (3) 

where the gradient V is with respect to the x variable. 

The corresponding dimensionless random lensing map r\ g : L — > R 2 , (the lens plane), is defined by (see 

r) g (x) =x- Vip g (x). 

The advantage of the lensing map is that solutions x of ^ are equivalent to solutions x of the lens equation (see 
p. 81): 

Vg(x) = y, 

which are preimages x of y under the lensing map n g . Solutions x € L of the lens equation are called lensed images. 
Microlensing also causes light rays to bend. For a light ray passing through the point x in the lens plane, the 
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bending angle a g (x) at x is given by (see 



21 



, p. 79): 



a g (x) = Vip g (x). 



Assumptions and Notation 



We shall abide by the following throughout the paper: 

• Rectangular coordinates (u, v) are assumed in the lens plane L. 

• All masses are equal, rrij = m =/= 0, for j = 1, • • • , g. 

• The star positions £* = (Ui,Vi) are independent and uniformly distributed in the disk B{0,R) of radius R 
centered at 0, that is, 

(Di,V 4 )~Unif {(«,«)€£: |u| 2 + \v\ 2 < R 2 }, i = !,-■■ ,g. 



• The quantities R and g are related by the following physical formula for the surface mass density k* : 

m g 

«. = -R2-- 

Unless otherwise stated assume that k* is fixed. 

• Notation: For clarity, the statement 

"-F(ar) = G{x) + 0(AT(a;)) as x -* oo" 



means precisely that 



lim l^)- G (,)| <c 

for some constant C > 0. 
• Notation: log 1 A = (log A)*. 

III. RESULTS ON THE RANDOM TIME DELAY FUNCTION 

We can write the time delay function as 

T g , y {x) = \\x-y\ 2 - ^M 2 + \{x\-x\) + J2W gJ (x), 
where 

w gA x ) =-mlog|a;-^| 
is the point mass potential. 



We define a "normalization" T* of T SiI/ as follows: 



Tl y ( X )^T g , y ( X )+gm\ogR=±\x-y\ 2 -^-\ X \ 2 + l(xl-xl) + J2 W IA 



where 






W* -(a;) = —to log 



(4) 



is the normalized point mass potential. Note that this normalization (of the time delay function) does not change the 
lensing map and relevant physical lensing quantities such as magnification, difference in arrival times between two 
images, etc. 



A. Exact Moments and Probability Density Functions 



The first and second moments of a random point mass potential are given by: 

Proposition 1. Let x <G B(0, R) and y el 2 (the light source plane). Then 

PIW , „ n.4R-\x\ni-21og(R-\x\)) 2k, f R+lxl .. . .... . , 

E Wg,j{x)\ = / rf(r, \x\)(logr)dr, 

2 9 Kg J R _\ X \ 



EK{X)] = mn.(R - \x\? [1 + 2(log(fl - N) - 1) log(* - |x|)] + 2^ r^N 



where f(r,\x\) — Arccosf ' x ' 2 , r , r — J for non-zero x, and f(r, 0) = 0. 



Proof: See Appendix. □ 

Both integrals in the above proposition arc finite (and can be computed numerically). Hence, the expectation of 



l\ X -y\*-%\ X \* + l(~* ^ 



E[T g , y (x)] = -\x y\ 2 ^\x\ 2 + Ux\ xj) + gE[W gj (x)} 



and its variance 



Var[T s ,,(z)] = gE\]V^(x)] - g(E[W gJ (x)}) 2 

are also finite. Instead of pursuing exact expressions for the mean and the variance of T* (x) for general iGfl(0,iJ), 
we take an asymptotic approach: 

Corollary 2. Let x G -6(0, R) and y GM. 2 (the light source plane). In the large g limit, we have 

E[T*Jx)] = ^g + 0( 3 1/2 log 5 ) 



Var[Tl y (x)} = (^y 9 + 0(g^ 2 log 2 g). 



Proof: See Appendix. 
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Corollary [5] shows an interesting link that has a "stochastic Morse theoretic" flavor, namely, at leading order, 

777 / 777 \ 

E[T;jx)} = -Si (A,), VarlT^^)] = (-J B X (X>*) 

where D^ — B (0, i?) — {^ 1; • • • , £ g } and Bi(D^) is the first Betti number of £>* . This is because Bi(D. t ) = g, which 
is the number of "holes" in D* (e.g., [21j, p. 399]). 

We now derive the p.d.f.'s of the point mass potentials W g j(x) and W* Ax): 

Proposition 3. For x G B(0,R), the p.d.f. ofW g j(x) is given by: 



fw gd (x){h) = < 



mW eX Pl mJ' 

A R Jh), 
0. 



- mlog(i?- \x\) < h 

mlog(R+ \x\) < h < -~m\og(R- \x\) 

h < -m\og(R+ \x\), 



where 



1 - [r h cos" 1 9 RA {h) + R 2 cos" 1 9 Rft (h) - B R , x (h)j /(irR 2 ), 

exp [2/7/777], 

(\x\ 2 +rl-R 2 )/(2\x\r h ), 

(\x\ 2 +R 2 -rl)/(2\x\R), 



A R , x (h) 

e R ,x{h) 

R , 2 (h) 
B R , x {h) = y/(R + r h - \x\){d + Th -R)(d + R- r h )(d + r h + R)/2 



(5) 



It then follows that the p.d.f. ofW*Ax) is: 



fw* d (x)(h) 



-exp[-% 



A lf (/7), 

0, 



-mlog(l-J^) <h 

777log(l + M) < h < -mlog(l - i^.) 

h< -mlog(l + J^). 



(6) 



Proof: Sec Appendix. 



□ 



At the origin, we obtain the exact p.d.f. for the random time delay function: 
Corollary 4. The random variable T* (0) has a shifted gamma distribution with density 



( 2 \ 9 {h-cY- 1 

\ m ) (g-l)\ 



exp 



Mh'o) h>c 

0, h < c, 



where c = \y\ 2 /2. 



Proof: From Proposition [3l we know that W„,-(0) is an exponential random variable with parameter 2/m. By inde- 
pendence, £?=i W^.(O) ~ r(«?, 2/m). D 

A direct calculation using the result of Corollary |4] shows that the exact expectation and variance of the normalized 
random time delay function at the origin are given respectively as follows: 



m _. rm „ ,_., (m 



2 



E{Tl y (0)} = - g + c, Var[T 9 yo)] = [-) g, 

which is consistent with Corollary 

Although Corollary 0] only gives the p.d.f. of T* (0), it will play a central role in establishing the asymptotic p.d.f. 
of T*[x) for an arbitrary x G B(0, R) — see Theorem [5] below. 

B. Asymptotic P.D.F. for the Time Delay Function 

The expression of the p.d.f. fx* t x ) f° r arbitrary x can be derived via convolution and equation ©, but is very 
cumbersome. As g — ► oo though, the analysis simplifies. 

Theorem 5. For every x £ B(0,R), let Jt* ( x ) oe ^ e P-d.f. ofT*(x). In the large g limit, we have: 



fT Sy {x){h) = f Tg%{0) (h-d 1 (x,y)) + O(g-^) 



9iV 



(A) 



(h-di-cY 

— - — exp 



2(h-d 1 -c) 



(9-1) 

0, h < d\ + c 



h> di + c 

0(g~ 3/2 ) (8) 



where d\{x,y) = \\x — y\ 2 - ^f\x\ 2 + %(x 2 - x\). 
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Proof: See Appendix. □ 

Theorem [S] shows that, at leading order, the p.d.f. fa* O) of the normalized time delay function T*{x) is a gamma 
distribution. Moreover, since 



h — d\ [x, y) — c = y —m log 



I* - 



R 

the p.d.f. fx* ( x ) for any x G B(0, R) depends at leading order on the normalized potential time delay due purely to 
point masses. 

IV. RESULTS ON THE RANDOM LENSING MAP 

The lensing map r\ g : L — ► M 2 is defined by 

r) g (x) =x-Vtp g (x), 

where L = R 2 — {£1, . . . , £ ff } is the lens plane. The components of the lensing map are given by 

, \ ,-, ^ v^ tnCUj —xi) ,-,,-, n v^ w(V" 7 - - x 2 ) 

m,g( x ) = (1 - «c + l)xi + }_^ H273 — > V»,gW = (1 - *c - 7)^2 + ^ H273 — 

j= i ■"•jW j=1 fl jW 

where i? 2 (x) = (t/j - xi) 2 + (V- - x 2 ) 2 . 



A. Exact Expectation and Variance of the Lensing Map Components 

The components rji tg (x) and r)2, g {x) of the lensing map have the following basic statistics: 
Proposition 6. For x = (x\, x 2 ) G Lf] B(0, R), the exact expectations of the components of the lensing map are 

E [Vi,g(x)] = [1 - («c +«») + 7] *i and .£[772,9(0;)] = [1 - («c + «*) - 7] *2. 

Boi/j variables have infinite variance. 



Proof: See Appendix. □ 
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Proposition [6] provides statistical information about the position values of the lensing map in the light source plane. 
In fact, suppose that an image position x = (0:1,2:2) is given. The proposition gives: 

E[r)i, g (xi)] = [1 - (k c + k*) + 7] xi, E[r)i >g (x2)] = [1 - (« c + «*) - 7] x 2- 

If x = 0, then on average the values of r] g (0) are isotropically distributed around the origin in the light source plane. 
More generally, we obtain the following: 

• Case I ( "macro-minimum" ) : 1 — (k c + k* ) + 7 > and 1 — (n c + k*) — 7 > 0. If a; lies in a given quadrant of the 
lens plane, then on average the random lensing map will send x over to the same quadrant in the light source 
plane. 

• Case II ("macro-saddle"): 1 — (k c + k*) + 7 > and 1 — (k c + k*) — 7 < 0. If x lies in quadrant I of the 
lens plane, then on average the random lensing map will send x over to quadrant IV in the light source plane. 
Similarly, we obtain results for the other quadrants. 

• Case III ( "macro- maximum" ): 1 — (k c + r* ) + 7 < and 1 — (ft c + ft* ) — 7 < 0. If 2: lies in given quadrant of the 
lens plane, then on average the random lensing map will send x over in the light source plane to the quadrant 
symmetric with respect to the origin (e..g, quadrant I goes to III on average). 

B. Asymptotic P.D.F. of the Lensing Map 

For the remainder of the paper, we shall assume that 

R 2 =9/k, 9>1, 

which gives ft* = nm and ignores the known case with g — 1 (single star lens). Also, we suppose that the random 
lensing map has an absolutely continuous cumulative distribution function. Then: 

Theorem 7. Let x = (.xi,:^) £ LP\B(0,R) be fixed. The p.d.f. of the scaled lensing map r '° is given in the large 
g limit by: 



(i.-ii)^-n-i,i- 



/%w(M) = 



/TSFs (y2-Ku g 



x\{h - a\) + x 2 (k — 5 2 ) 



\/log5 o-g 



k\ {{h-~ ai ) 2 + {k-~a 2 f)-2dl log(log 5 ) 
4 



4?r <5-4 log g 



O' (9) 
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where (h, k) are the possible values of 'l? and 

(1-k c + 7)xi . (l-«c-7)*2 



«i 



Vlogg 



«2 



Vlog.9 



Vlog5 



T g = —^= ylog(B g 1 / 2 ), B = (7,, is i/ie Euler constant). 



>n 



Av* 



The p.d.f.s of the components of the scaled lensing map are: 

(g-M 2 



fii, g <.*){£) 



'2n a„ 



1- 



k* aii(c-Ot) k 2 (c-ai) 2 -^ 2 log(logp) 



Vlog.9 ol 



4tt 



3 



log 5 



O 



log.' 



for i = 1,2. For notational simplicity, we use c to represent the possible values of both rji >g (x) and r)2, g {x). 



Proof: See Appendix. 



D 



Observe that the p.d.f. of the scaled lensing map f ■ngW depends on x in its leading order term. Furthermore, the 



fact that 



Var 



m, g ( x ) 



00, 



« = 1,2, 



Vl°g9_ 
makes applying a standard Central Limit Theorem unwieldy, thus the method of proof presented in the Appendix. 

Theorem [7] shows that at leading order, the p.d.f. of the scaled lensing map is bivariate normal and illustrates how 
at the next order the p.d.f. deviates from normality. Moreover, since 

(h -ai,k- 62) = "i— - = a* 
Vlogg 

which is the scaled bending angle due purely to stars, the quantity (h — ai) 2 + (k — c^) 2 in (|9|) is 



U,*|2 



In other 



words, the leading factor of the p.d.f. of the scaled lensing map depends on the square of the magnification of the 
scale bending angle due to stars. This radial symmetry is broken at the next order term since there is a dependence 
on the individual components of the scaled bending angle. 
We now discuss several consequences of Theorem [7J 



Corollary 8. As g — > 00 with fixed iel ! , we have: 



( 



Vhg( x ) m,g( x ) 

Vlogg' Vlog.9 



(Vl,oo,T]2,oo)- 
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Here "=$■" denotes distributional convergence and (?7i,oo, ?72.oo) is a bivariant normal random vector with independent 
mean-zero normal random variables as components, each with variance nl/(2ir). 



Proof: See Appendix. □ 

We can readily obtain the p.d.f. of the unsealed lensing map: 
Corollary 9. Let x — {x\,X2J € L be fixed. The p.d.f. ofr) g (x) is given in the large g limit by: 



ih-a i y + (k-a 2 ) 

e ^ 

fri„(x)(h,k) = 



x ± (h - ax) + x 2 (k - a 2 ) 



-I 



K 2 ((h-a 1 ) 2 + (k-a 2 ) 2 -2o-2) 

- I I \ I t i V / \ 

log 9 



Air a* 



-B{g) 



^(rW (io) 



with the possible values of the random vector T] g (x) written as [h, k). Here a g and H(g) are as defined in Theorem^ 
and 

ai = (1- k c + j)xi, a 2 = (1 - k c - 7)«2, B{g) = log(logg). 



Proof: First, note that the possible values of r) g {x) and SpS are related by 



h = — . k 



yfiogg Vlog 

We know that 

1 . ,r rs 1 . ( h k 



1 - - 1 / 

f Vg (x)(h,k) = / „ aW (h,k) = /„„(,) - 



gff Vi5g» Wlog5 vlogl 
The result then follows directly from equation ^j. □ 

Let us consider the leading function in equation (|10|) . namely, 



14 



(h-a 1 )' 1 + (k-a 2 ^ 

^71 



fg,x(h,k) 



2-ko, 



x 1 (h~a 1 )+x 2 (k~a 2 ) nl (g» ~ a ^ + ( k ~ a *? ~ 2 °t) 



where x £ B(0,R) and (h,k) £ M. 2 . Note that H(g) > if g > 3. The function T g , x is actually a p.d.f. for 
\x\ sufficiently small and g > 3. In fact, under the latter, we have J r 9}X (h,k) > for all (h,k) £ M. 2 . Moreover, 
integration-by-parts gives: 



/ J-g. x (h,k) dhdk = 1. 

Js. 2 



Hence, the function T g>x is a p.d.f. on M 2 for \x\ sufficiently small and g > 3. Figure Q] depicts the graph oiT g ^ x with 
x = (x 1 ,x 2 ) = (0.2,0), k c = 0.405, 7 = 0.3, k„ =0.045, and fif = 10 6 . 




FIG. 1: Asymptotic p.d.f. of the random lensing map in the large number of stars limit. 

Furthermore, there is enough control in the p.d.f. f n ( x ) to perform the expectation integrals to obtain: 

Corollary 10. Let x £ L the lens plane. In limit g — > oo, the expectations of the components i]i,g(x),r]2,g(x) of the 
lensing map are given by 

E[rn, g (x)] ~ [1 - (k c + k») + (-l) l+1 7 ] Xi, i = 1,2. 

j4Zso ; i/ie variances of the lensing map components are infinite: 



Var[^ s (x)] ~ -lpog(B 5 1/2 ) + - log log 5 - tt*?] = a* 



--log log 

Z7T 



2 2 
K * x i, 9 ~~ * °°i 



i = 1,2. 
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Proof: This result is obtained from Corollary [9] (see equation (flC)]) ) by computing the expected value and the variance 



of the random vector with density F g>x (h, k). 



□ 



Corollary 1101 shows that the expectation from the asymptotic p.d.f. T g _ x is consistent with the findings in Propo- 
sition [6) We also have consistency for the variance as g — * oo. Note that in the leading order of Corollary [TOl the 
divergent variance depends on the position x. 

Next, observe that given g > 1, the normal leading factor of the p.d.f. /„ r x \ is maximized at 

a = (ai, 02) = ((1 - k c + 7)11, (1 - « c - 7)2:2). 

The following explores the probability that rj g {x) is located in a neighborhood of the maximum: 

Corollary 11. Let x — (x\,X2) € L, fix ro € R, and let n be a positive integer. The probability that the distance 
between r) g (x) and the point a = (ai, 02) is between (n — l)ro and nro, is given in the large g limit by 



P 



(n - l)r < \Vg(x) - a| < nr 



eM-^} 



4ttct4 



e"» (47r^ + (n-l) 2 r^^loglog 5 ) 

- e^(4iro- 4 g + n 2 rlnl loglogg) 



O 



log 9 



11) 



Proof: See Appendix. 



□ 



Notice that the displayed leading terms in equation (jlip are independent of the position x, continuous matter 
density k c , and external shear 7. 

Corollarv llll combined with the lens equation, provides probabilistic information about the angular position of the 
light source. Indeed, for n = 1 the corollary gives the probability of T) g (x) lying in a disk of angular radius ro centered 
at a, while for n > 1 we obtain the probability of Tj g (x) lying in an annular region between radii (n — l)ro and nro 
centered at a. For illustration, we compute these probabilities in some specific cases: 

Po.i(n) = P[(n- 1)/10 < |»fc(aO - a| < n/10 

and 



Po.i(n) = P[0 < !%(») - a| < n/10] = ^Po.iW 
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for n = 1,2,3 with the constants rg = 0.1, k* = 0.045, and g — 10 6 . The probabilities are given in Table I. 



n 


1 


2 


3 


Po.i{n) 


0.5595 


0.4063 


0.0337 


Po.i{n) 


0.5595 


0.9658 


0.9995 



Table I: Probabilities of a lensing map's value lying in various annuli (second row) and disks (third row). 

The second row of Table I shows that the lensing map value n g (x) has a 56% probability of lying inside the disk of 
angular radius ro = 0.1 centered at a, but a 97% probability of lying within radius 2rQ. The probability of being in 
the outermost annulus 2ro < \r) g (x) — a| < 3ro drops dramatically to 3%. 

V. RESULTS OF THE RANDOM BENDING ANGLE AND FIXED POINTS 

We now use our study of the lensing map to obtain probabilistic information about the bending angle vector and 
lensing fixed points. 

A. Asymptotic P.D.F. of the Bending Angle 

Recall that the bending angle vector a g {x) is given by a g (x) — '\7ijj g (x). Its asymptotic behavior is given below: 
Corollary 12. Let x — {x\, X2) G L. The p.d.f. of the bending angle vector is given in the large g limit by 



^ 



fa g (x){h*,k*) 



(v^<? g y 



1 — K 



xi(h* — a|) + X2(k* — a^) 



c 53 log(log9) ] + °(^- (12) 



where (h*,k*) denotes the possible values of the random vector a g (x) and (a*,^) = (( K c — 7)^1, ific + 7)^2, while 
o~ g is defined in Theorem [^ In other words, for sufficiently large g, the bending angle is approximately a bivariate 
normal whose components a\ tg {x) and a2,g(%) are normal random variables with respective means 



[—(k c + «*) +7] x\ and [-(k c + ««) -7] »2) 
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and variances diverging at leading order as 

n 2 i \ n 2 1 

-^[log(B 5 1/2 ) + -loglog ff -7ra;f] and -^[\og{B g 1 ' 2 ) + - log log ff - ttx 2 2 ], 

■K I 7T 2 

respectively. 

Proof: By definition, we have a g (x) = x — T) g (x) — x + N g (x), where 5; = ((— k c + j)xi, (— k c — 7)3:2). The corollary 
then follows from properties of independent random variables, Theorem [71 and some of its corollaries. □ 

When x = 0, the second part of Corollary [12] recovers, in leading order, earlier results by Katz, Balbus, and 



Paczynski [7j and Schneider, Ehlers, and Falco (see [23J], p. 325). Also note that from the discussion after Theorem[7J 
we have that the leading order term depends only on the squared magnitude of the bending angle due purely to point 
masses, and the next term breaks this radial symmetry. 

B. Fixed Points 

Lensed images of a light source usually do not appear at the same angular position of the light source. In the 
case where they do, they are called fixed points, which is a notion introduced in gravitational lensing by Petters and 
Wicklin 20( . More precisely, fixed points of the lensing map r\ g are points igL such that 

n g (x) = x 

or equivalently, points where the bending angle vector vanishes. a g (x) = 0. Since a g is a continuous random vector, 
we have P[a g (x) = 0] = 0. Nevertheless, we can study how close n g (x) is to x. 

Corollary 13. Let x G -B(0, R) C L and e > 0. The probability that rj g (x) is within e from x in the light source plane 
is given by 

P[\n g {x)-x\<e]= f f aq{x) (h\k*)dh*dk*. 

J B(0,e) 

Proof: This follows directly from Corollary [T2l and the above discussion. □ 
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VI. CONCLUSION 

We presented first steps in the development of a mathematical theory of stochastic microlensing that focused on 
the building blocks of the theory, namely, the random time delay function and random lensing map. We derived exact 
analytical formulas for the expectation and variance of the random time delay function and random lensing map about 
any point. In the large limit of stars, we found a simple asymptotic expression for the expectation and variance of a 
normalized time delay function at an arbitrary point. For the same limit, we highlighted an interesting link between 
the leading order term of the expectation and variance of our normalized random time delay function and the first 
Betti number of its domain. In addition, the asymptotic p.d.f.s of both the normalized time delay function and the 
scaled lensing map were characterized in the large number of stars limit: The asymptotic p.d.f. of the former is a 
shifted gamma density at leading order, while for the latter the p.d.f. is a bivariate Gaussian distribution. The p.d.f. 
of the random scaled lensing map is also shown to depend on the magnitude of the scaled bending angle due purely to 
point masses at leading order and we illustrated explicitly how this radial symmetry fails at the next order. We also 
derived the asymptotic p.d.f. of the random bending angle vector and gave an estimate of the probability of a lens 
plane point being close to a fixed point. Overall, the paper determined and illustrated new analytical results about 
the microlensing behavior of random time delay functions and random lensing maps about an arbitrary point of the 
lens plane. 

In Paper II of this series, we shall explore the microlensing random shear and expected number of microimages. 
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APPENDIX A: PROOFS— THE RANDOM TIME DELAY FUNCTION 

Proposition 1. Let x e B(0, R) and j/£R 2 (the light source plane). Then 

W\W ( M ^-H) 2 (-l + 21og(i?-H)) 2«. [ R+M 

E Wg,o\ x )\ = o / r -f( r ' m){logr)dr, 
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and 



E[WlAx)\ = 



mK*(R - \x\) 2 [1 + 2(log(i? - |x|) - 1) log(R - |i 



2ff 

2m K * f R+ W 



n 9 Jr-\x\ 
where f(r, \x\) — Arccosl ' x ' 2 m^~ — ) / or non - zer o x, and f(r, 0) = 0. 



r / ( r > I ^ |) (log 2 r)dr, 



Proof of Proposition 1: Let x — (xi,X2) € B(0,R) n L and y e R 2 . Translate the rectangular coordinates (u, w) so 
that its origin moves to position x. Denote the resulting new coordinates by (u',v') and its origin by . In the 
system (u',v'), the old origin now has coordinates (u',v') = —X. Define luq to be the unique principal angle with 
cos wo = xi/|x|, sinwo = X2/\x\, and let uj = luq + tt. Rotate the rectangular coordinates (u',v') counterclockwise 
by angle u> to obtain new coordinates (u",v") with origin 0" = 0'. Note that the old origin now lies at position 
(u",v") = (|x|,0) = x" on the positive w"-axis. Finally, let (9,r) denote polar coordinates in the frame (u",v"). The 
figure below illustrate some of these change of variables. 





\ 










Ll. 1 j^ 




s\ x 








n 




I u 


u"V/ 








x" 



FIG. 2: Change of variables. 



We have 



E[W gd {x)] = -mE[log\x-Z J 



20 



and 



E[log\x-^\] = -^ f log[( 

^ti Jb(o,r) 



u — X\) + (v — X2) ]dudv 



2irR 2 



B(o',R-\x\) 



log [{u'f + {v'f} du'dv' 



I B(-x,R)\B(o' ,R-\x\) 

For x = 0, the result follows. Suppose i^O. Then 



log [{u'f + {v'f} du'dv' 



E[\og\x-^\] = 



B(0",R-\x\) 



2irR 2 



IB(x",R)\B(0",R-\x\) 
1 r R-\x\ , r-2iT 



log [{u"f + (v"f] du"dv 



*R 2 J a 

2 rR-\x\ 



(log r)( de)dr+- 



log [(w") 2 + (v") 2 ] du"dV 
I r R +\ x \ 



R 2 



R-\x\ 



f(r,\x\) 

r(logr)( / d6)dr 



-f(r,\x\) 



R 2 J a 



2 rR+M 
r log rdr -\ — j rf (r,\x\) log rdr 



irR 2 



R-\x\ 



(R-\x\) 2 (-l + 2log(R-\x\)) 2 f R+ ^ 



2R 2 



which gives the mean. The second moment is obtained similarly. 



ttR 2 



rf(r,\x\)(logr)dr, 



R-\x\ 



D 



Corollary 2. Let x G B(0, R) and y £ M 2 (the light source plane). In the large g limit, we have 



E{Tl v {x)] = 1 ^g + 0{g 1 / 2 \ogg) 



and 



V & T[Tl y {x)] = ^g + 0{g 1 ' 2 log 2 g). 



Proof of Corollary 2: Let d n (x,y) — (^\x — y\ 2 - ^-\x\ 2 + ^(x 2 — x^U/n, n e N for arbitrary x. 



E[T;jx)] = d 1 (x,y)+gE[W giV (x)}+mg\ogR 

HtI{2[ 2 -m^v) + (i-\x\/R) 2 -2- 



2|x|i?logi? 



K*R 2 

-2(l~\x\/R) 2 \og(l-\x\/R) 



R R 

4 fR+\x\ 



ttR 2 



R-\x\ 



rf(r,\x\)(logr)dr 
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But 



R+\x\ 



R-\x\ 



rf(r,\x\) (log r)dr 



rR+\x\ 

< / \rf(r, |x|)logr| dr 

JR-\x\ 

<-R+\x\ 



/■/t-t-|X| 

< 2ir / r (log r)rfr for sufficiently large R 

JR-\x\ 



>R-\x\ 

27ri?logi? 



so 



-R+M 



R-|x| 



r/(r, |x|)(logr)dr - O(RlogR) = O^ 1 ' 2 log g). 



The first result then follows using the relation k* = mg/R 2 . For the second: 



Var[T* (x)] = gVax\W g , v {x)] 



= gE[Wl y (x)}-g(E[W g , y (x)}f 

(R - \x\) 2 [l + 2(log(i? - \x\) - 1) log(R - \x\ 



(R-\x\) 2 (-l + 2log(R-\x\)) 



0{R\og z R) 

2 



O(RlogR) 



mK*(R — \x\Y 



i (R-W) 2 

2 4R 2 



+ (log 2 (R-\x\)-log(R-\x\))(l 
V^^ + (R\og 2 R). 



(R-\Af 

R? 



0{R\og 2 R) 



This completes the proof. 



□ 



Proposition 3. For x <G B(0,R), the p.d.f. ofW g j(x) is given by 



fw g ,j(x){h) 



2 r_2fci 

mix r L m J ' 



A R Jh), 

0, 



- mlog(i? - \x\) < h 

— mlog(R + \x\) < h < —m\og(R — \x\) 

h < -mlog(R+ \x\), 



where 
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A R , x {h) = l-(r h cos- l e R ^{h)+R 2 C o^ 1 9 R . 2 {h)-B Rx {h))/{nR 2 ) 1 

r/, = exp[2h/m], 

6 R1 (h) = (\x\ 2 +r 2 -R 2 )/(2\x\r h ), 

0R,2(h) = (\x\ 2 +R 2 -r 2 h )/(2\x\R), 

B R , x {h) = ^/{R + r h - \x\){d + r h - R){d + R- r h )(d + r h + R)/2. 



It then follows that the p.d.f. ofW*Ax) is: 



l\Y 



;» 



(h) 



-exp[-^], 



A hj _(h), 
0, 



-mlog(l-^) < h 

mlog(l + J^) < h< -mlog(l - §) 

ft< _ m iog(l + H). 



Proof of Proposition 3: Given x £ B(0, R), consider the point mass potential 



w g,j( x ) = -mlog|a;-|j|. 



Then 



P(W g ,j(x) <h) = 1 - P[\x -£\< exp(-tym)] 
axe&(B(0,R)nB(x,R)) 



1- 



irR 2 

cxp( — 2h/m) 



(see 



A Ri \ x \{h), 
0, 



< exp(-/i/m) < R - \x\ 

R—\x\< exp(-/i/m) < R + \x\ 

R+ \x\ < exp(— h/m), 



(Al) 



(A2) 



22] for the derivation of the middle entry). Thus, the function P(W g j(x) < h) is piecewise smooth in h, and 



therefore, for almost all h we obtain: 



d_ 

dh' 



fw gA *){h) = it p lw g Ax) < h], 



which yields (0 and ©. 



□ 



2: J > 



Theorem 5. For every x G B(0,R) 7 let f T * (*) &e i/ie p.d./. ojT* {x)/g. In the large g limit, we h 



/*■;,„<.) (ft) = /- 9 «, y <o) (ft - d„(x, y)) + 0(.g- 1/2 ). 



Proof of Theorem 5: Let /ii = — mlog(l + \x\/R) 1 /12 = — mlog(l — \x\/ R), and <pT* ( X ) be the characteristic function 
oi Ilv^l_ It follows that: 



fe-^^'V^wWl =-/ e ah ^ e - 2h / m dh+ e lth/g f w * {x ){h)dh 

= (i _ ^-1 _ 1 [ h2 eW g -^)h dh 



+ / e utl ' g f wl . (x) {h)dh 



In 



^, 0) {t)+0{g-^l 



where we used the following facts: 



hi 



2/f 



fw* ti (x){h)dh 



P[W* gJ (x)<h 2 ] = l-(l-^) 2 , 



Z-% h dh = (l-^) 2 -l. 

R 



mg J 



he~™ h dh = - 



it 



9 Jh! 



h 2 (i - p[w:ax) < h 2 \) - -p[w:ax) < h 2 



hfw* .( X) {h)dh = - \h 2 P[W^(x) < h 2 \ + 0{g- 1 ' 2 ) 

9,3' ri a '•* 



h 2 --P[Wl {x)<h 2 ] = Oig- 1 ). 



Therefore, 



e- M ^y^ Tlv(x) (t) = ^ TJ ^{t) l + 0(.g- 3 / 2 ) 

a 9 

= y^ y( 0)(t)[l + O(g' 1/2 ) 

where the expansion series (sums of terms of order less than g 1 ' 2 ) is the sums of terms of the form g q x t n with 
q e Q n (-00, -1/2) and n e N. 
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By Proposition [3l the distribution function of T* (x) is absolutely continuous. Therefore, using Inverse Fourier 
transform, we obtain 



M. H (*> W = /^. H c°) (h - d g (x, yj) + Q(g- 1/2 ). (A3) 

9 3 

We can now use the following fact to complete the proof: 

where c > and X is a real- valued random variable. □ 

APPENDIX B: PROOFS — THE RANDOM LENSING MAP 

Proposition 6. For x = (x\,X2) € LC\ B(0,R), the expectations of the components of the tensing map are 

E [Vi,g(%)] = [1 - (« c +«*) +7] *i and 2?[?72, 9 (x)] = [1 - (« c + «*) - 7] X 2 . 

i?ot/i variables have infinite variance. 

Proof of Proposition 6: Let x — {x\,Xi) G B(0,R) ("I £ and y <E R 2 . Let /(r, |x|) be defined as in the statement of 
Proposition Q] and employ the coordinates (u',v') and (u",v") defined at the start of the proof of that proposition. 



gm f (u — Xi) 

kR 2 Jb(o,r) ( u ~ x i) 2 + (v- x 2 ) 2 



E[Vi, 9 (x)} = (1 " «c + 7)^1 + -pa / , , 7-— r^— 7T— TT2 dudw 



(1 - k c + 7)^1 + — / nwn A2 du ' dw ' 

77 7s(0,ii)\B(0',fl-|x|) K) +(«) 



If a; = 0, the result follows. For x 7^ 0, we have: 



," ^,^0 , . 1 „," , 



/" u' , , , , f W"COSW + U"SlHW ,//,// 

/ , f 7 ,A2 , r „A2 du dw = / (,.„S2 f.,,,\2 du dv 

J B(0,R)\B(0' ,R-\x\) V u ) + ( v ) Jb(0,R)\B(0",R-\x\) ( u ) ~ K v ) 

rR+\x\ rf(r,\x\) 
= I I (cos 9 cos lo — sin 9 sin uj) d6 dr 

jR-\x\ J-f{r,\x\) 
pR+\x\ 

— 2(cosw) / sm[f(r, \x\)]dr 

JR-\x\ 



t R+ \ x \ / /Id 2 + r 2 - i? 2 \ 2 

2(cosc) / wi- r 1 1, 1 ' dr 
Jr-\x\ v v 2r M 

7r|x| cosw 



= —-KX\. 
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Therefore, E[rji^ g {x)\ = (1 — n c — k* + 7)0:1. Similarly, E[r]2,g{x)] = (1 — k c — k* — 7)0:2- 

Next, we establish only the divergence of the variance of rj\^ g {x) since the proof for the case of r/2,g(x) is similar: 



Var[?7i,g(x)] = ^Var[ 



m(Uj — xi) 



3 = 1 



(u J - Xl r + {v -x 2 f 



> 



gm 2 f u 2 gm 2 2 



ran* II sm " ,„ mK^x 

hm / r / — - — dOdr 



IT e^O 
= 00. 



R 2 



D 



Theorem 7. Let x = (x\, x 2 ) E L be fixed. The p.d.f. of the scaled lensing map v y x > is given in the large g limit by 



logs 



f ng(x) (h,k) = 


(V^z 2 ) 

4ir 


k* xi(h ~ di) + x 2 (k - a 2 ) 




Vlog5 CTg 

((h-d 1 ) 2 + (k-d 2 ) 2 )-2a 2 log(log 5 )" 

^ lOg ff 


+ (l' 



with the constants as previously defined. The p.d.f.s of the components of the lensing map are: 






/^jW (5) = 



'27T(T 



1- 



K* Xj(£ — Oj) K^ (5 " 


-a,) 2 - 


-a 2 log(logg)- 


v log5 


Vlog.9 cr2 4tt 


^ 


log5 - 



fori = 1,2. 



Proof of Theorem 7: The scaled lensing map can be written as: 

r)g(x) __ N 9 (x) / (1 - K c + 7)xi (1-K C - 7)0:2 



Vlog.9 V lo g.9 



Vlog5 



Vlogg 



where 



Vlogg 



«i, a (a;) N 2 , g (a:)) 



/logs ' VTogg 



with 



(Bl) 



Ni, ff (*) = £ 



m(C/j — xi) 



=1 W 



N 2lS (x) = ^ 



m(V^ - 0:2) 



=i w 



The possible values of the random vector ^pSj will be written (/i, fc), while those of ,y- x ' will be denoted by (h, k). 
We begin by determining the asymptotic p.d.f. of ^ x ' . Consider an arbitrary x — (xi,x 2 ) e B(0,R). For 
the random vectors (Ui,Vi), set u x = u — 0:1, v^ = f — 0:2, and t 2 , = u 2 ^ + v 2 2 . The joint characteristic function 
y »,w fti)*2) of 4j^=i satisfies 



,. . C\ l f r im f tmx + t2Vx\-\ 

y %M (tiM) = —52/ exp ^ 

2 [ R ~ lxl rJ /mVtT+tf,-!- 



dudv 



R 2 Jo 



where Jq is the zeroth Bessel function. 



Vlog.g 



') 



+ 



nR 2 



exp 



S(0,ii)\S(a;,fl-|a:|) 



I" 177 

Ufa. 



im ftiUx +hv x 



V^ogg 



dudv, 



Adding and subtracting -j^ J Q rdr to the first term in (|B2|) yields: 



2(i 



(B2) 



2 rR-\x\ 
Wo 



rJo ( Wg+V iy,,. 



viog 



= l - 



2N 
R 

2\x[ 
R 



R-\x\ 



\x\ 2 2_ 

~W~R? Jt] 

R 2 4tt 

+ 0(. 9 - 1 (log. 9 )- 1 ), 



1-J< 



2 K 2 (t 1 +il)log( J B 2 5 )+loglog(? 



»(' 



?VlogS 



dr 



Slogs 



(B3) 



where -B = e , t- and we used k* = Trm. 
For the second integral in ()B2|) . we get: 



irR 2 



exp 



S(0,fl)\B(x,il-|x|) 

1 



Z771 ( t\U x + i2^£C 



Vlog'3 



dudv 



nR 2 



I B(0,R)\B(x,R-\x\) 

2\x\ \x\ 2 i\x\n 



1 



im f t\u x + £2^3 



%/log.g 
^(i 1 ,t 2 ) + 0(. 9 - 3 / 2 (log. g )- 1 ) 



dwdw +0( £ T 3/2 (log 5 )- 1 ) 



R R 2 gVWd 

where E(t\, £2) = t\ cos uj + £2 sinw with to as defined in the proof Proposition [5] and k* = nm. 
Combining (jB3|) and (|B4|) . equation (|B2|) becomes, 



(B4) 



<P W g (x) (tl,f 2 ) 



1/9 



= 1- 



««(*? + *l)log(B 2 3) i|a;|K» 



4?r glogff .g%/logg 



E{h,t 2 ) 



Kl{t 2 +t 2 )log\ogg 
4tt .glog.g 



+ O^-^logg)- 1 ). 



Therefore, 



27 



¥> w fl (») (h,t 2 ) 



l ( K l(ti + tl)\og{B 2 g) 
4tt logg 



1 - - 
5 



l(K 2 M+i%) log{B*g) 

log 3 



4tt 



»M«* P /. ,n ft|(tf + t|)loglogg 

— -=^=b(t 1 ,t 2 ) -: ; 

5 V log 5 4tt fflog^ 

+ 0{g- 1 (logg)- 1 )] 9 
^^-Etti t 2 ) - K *^i +t 2) lo S lo g5 



exp 



(4 + 4)a* g 



2 log g 



l\X\K* 



E(h,t 2 ) 



+ 0((log 5 )- 1 )" 
«2(<i+^)log(logs) 



log g 



O 



log 5 



(B5) 



\/log5 4tt log g 

The term remainder term in equation (|B5|) is a sum of term of the form p(g) x q(ti,t 2 ) where p(g) has order 

1/logg or less and q(t±,t 2 ) is independent of g and integrable. We can, therefore, take the Inverse Fourier transform 

(I.F.T.) of (|B5| to obtain the joint p.d.f of -^1. Since the I.F.T. of the L.H.S. is well defined due to r) g (x) having 



an absolutely continuous cumulative distribution function, the I.F.T. of the R.H.S. of equation (|B5p carries through, 
and can be computed term by term. The I.F.T. of the first term in (|B5|) is a standard integral, and gives 



1 



exp 



(2tt) 2 
For the second term in (IB5I). we obtain 



i(tih + hk) 



exp 



{t\+t 2 )<7 



2 log g 



h z +k z 



dt\dt 2 



2(^/lo gs ) 



2ir(a g /VWg)) 2 



(27T) 



exp 



-i{t\h + t 2 k) 



exp 



2 log .9 



IX /"£;); 
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where h cos v + k sin v = y /i 2 + fc 2 cos(^ + 0) , W. 
Similarly, 
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The sums gives 
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The I.F.T. of the third term in (jB5j) yields: 
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Finally, we combine these results to get: 
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Returning to J 3 ' ' in equation (|B1[) , we note that 
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The first part of the theorem now follows; the second part is a consequence of the first by integration. 
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Corollary 8. As g — > 00 wi£/j /raed ieM 2 , we have 

f m,g( x ) V2g(x) 
W lo g3' Vlog.9 



(^1,00,^2,00)- 



29 

Here "=$■" denotes distributional convergence and (r)i t00 ,T]2,oo) is bivariant normal random vector with independent 
mean-zero normal random variables as components, each with variance k 2 ./(2tt). 



Proof of Corollary 8: This result follows directly from equation (|B5[) in the proof of Theorem [7J 
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Corollary 11. Let x — (x\,X2) € L, fix vq S K, and let n be a positive integer. The probability that the distance 
between n g (x) and the point a = (ai, a?) is between (n — l)ro and nrg, is given in the large g limit by 
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Proof of Corollary 11: Set p ro ( n ) — P (n — l)ro < \Vg( x ) ~ a l < nr o 
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